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Abstract 

In the present paper we will introduce a new approach to multivariate 
interpolation by employing polyharmonic functions as interpolants, i.e. 
by solutions of higher order elliptic equations. We assume that the data 
arise from C°° or analytic functions in the ball Br. We prove two main 
results on the interpolation of C°° or analytic functions / in the ball Br 
by polyharmonic functions ft of a given order of polyharmonicity p. 
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1 Introduction and statement of results 

Interpolation theory is one of the oldest and most classical subjects of math- 
ematical analysis. It has been established in the work of Newton, Lagrange 
and numerous other mathematicians. Interpolation plays a fundamental role 
in algebraic geometry and numerical analysis, in particular in approximation of 
integrals (quadrature and cubature formulas), in finite clement methods, and 
others. 

There is a number of approaches to multivariate interpolation which are 
based on multivariate polynomials and radial basis functions (RBF), see e.g. 
[2D] , [5] , [TH] , [7] , [IS] . From the practical point of view the problem of interpo- 
lation of scattered data has been treated successfully by means of tools such as 
RBF (see e.g. [TH], [5] and references there) or polysplines (see [H]), which are 
in general not globally analytic. 

So far there remains the fundamental problem from the point of view of 
mathematical analysis to construct a multivariate interpolation theory based 
on globally analytic tools. The multivariate polynomials fail to deliver such 
tools. Indeed, it is well known and quite clear that multivariate polynomial 
interpolation differs in important ways from its univariate counterpart. The 
main difference is the fact that the multivariate polynomials fail to constitute a 
Chebyshev system, cf. [T], [TO] , 

Furthermore, let us recall that in the one-dimensional case the polynomial 
interpolation is closely related to a wide class of quadrature formulas. And the 
existing multivariate interpolation theories mentioned above do not provide a 
satisfactory theory of multivariate cubatures. 

On the other hand, objects like the solutions of elliptic PDEs, in particular 
the polyharmonic functions, have entered the scene of approximation and spline 
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theory (see e.g. [TT], [T2], [T3], [H], [5], [TS] and the references given there), and 
they satisfy a generalized definition of a Chebyshev system, see [16]. Is there 
an interpolation theory based on solutions of elliptic PDEs which provides a 
satisfactory analog to the classical one-dimensional results? 

In the present paper we address the above question by considering an inter- 
polation theory based on polyharmonic functions. Let us recall that a function 
h is polyharmonic of order p in a domain D C W 1 if it satisfies the equation 
A p h (x) = in D, cf. [2J, [2D]. It is important to emphasize the fact that in 
order to obtain satisfactory interpolation results one has to reconsider the whole 
paradigm of "set of interpolation points". In particular, in view of the fact that 
the space of polyharmonic functions is infinite-dimensional, one may consider 
interpolation sets T which are the union of hypersurfaces in 1™. Some results 
towards this interpolation theory have been obtained in [3], [H], [T2J, [8]. Let 
us focus on the analogy with the one-dimensional case: one is seeking such sets 
r which would correspond to the usual N points {xi,X2, ...,xjv} in K 1 where a 
polynomial P of degree < N — 1 solves the interpolation problem 

P{x i ) = c j for j = 1,2, ...,N 

for arbitrary data Cj. In particular, P (xj) = for j = 1, 2, N implies P = 0. 
It is clear that the main problem is to identify multivariate analogs to the 
"unisolvent" sets {xi, x 2 , x^}. 

Let us draw the reader's attention to the obstacles faced by the usual the- 
ory of interpolation with polyharmonic functions, related to the zero sets of 
polyharmonic functions. In [3], [12], [8], and references therein, such sets of 
interpolation T G D have been considered which are unions of N concentric 
spheres. It has been proved in these works that T is a set of uniqueness, i.e. if 
A N h — in D and if h (x) = for all x 6 V then h = in D. So far, attempts 
to consider sets F with a slightly more general geometry have led to a dead-end. 
In [3] (see the Russian edition of 1985) Atakhodzhaev has constructed a set of 
two closed convex curves 71 and 72 in IR 2 with contained in the convex hull 
of 72, such that there exists a (non-trivial) biharmonic function h with A 2 h = 
inside 72 and h (x) — for all x G 71 U 72 . This result has been dealt with in |8J 
as well. 

The last fact completely destroys any hope of finding reasonable unisolvent 
sets living in the space R™. In the present paper we formulate a concept of 
interpolation where the unisolvent sets live in what we call a " semi- frequency 
domain" which arises from the Laplace- Fourier spherical harmonic expansion of 
a function, see formula ([2J below. 

In order to motivate our approach to polyharmonic interpolation let us 
recall that in the classical one-dimensional interpolation theory error estimates 
are proved when data Cj are obtained from a differentiable function, i.e. 

cj = f{ Xj ) for j = l,2,...,iV, 

with / € C N+1 . In that case one may consider estimates of the error of inter- 
polation 

En [/] (x) = / (x) - P N (x) , 
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see [T7], [5]. More subtle results are obtained when / is an analytic function 
and N — ► oo. 

Now let us turn to the multivariate situation. Corresponding to the univari- 
ate case, in order to obtain a reasonable multivariate polyharmonic interpolation 
theory we will assume that the multivariate data arise from C°° or analytic func- 
tions. 

Let us first introduce some necessary notions and notations. We will work 
in the ball Br defined by 

Br := {x £ R" : \x\ < R} . 

Assume that we have a basis of the space of harmonic homogeneous polynomials 
of degree k (called spherical harmonics) which are denoted as Ykt (x) for k — 

0. 1, and £ = 1,2, dk, where 

d * = 1 \m (n + 2fc - 2) (n + fc - 3) • ■ • (k + 1) , (1) 
(n — 2)\ 

see [21]. They are assumed to be orthonormalized with respect to the scalar 
product 

— !— / u{6)v{6)d6 

on the unit sphere, where w„_i is the area of the unit sphere §™ _1 in R"; we 
have put 

x = r0, r — \x\ . 

Let us denote by C°° (Br) the set of C°° functions on a neighborhood of 
Br. For / G C°° {Br^J we have the expansion in spherical harmonics 

oo dk 

/(^EE/mM^mW- (2) 

fc=0 1=1 

We will use the following representation of C°° and of analytic functions in the 
ball, see [H p. 501, Proposition 1]: 

Proposition 1 Let f be in C°° (B_r) ■ Then we have the following expansion 

oo dk 

/(*)=EEMr 2 )r fe r M (0) (3) 

k=0 i=l 

where the functions f^e £ C°° ([0,i? 2 ]) . The function f is analytic in some 
neighborhood of in R™ if and only if there exist to > and M > such that, 
for all indices k > 0, 1 < £ < dk, and j > 0, we have 

< M k+j+1 j\. (4) 

1. e. f is analytic if and only if holds, and in that case each function fk.e is 
also analytic. 



sup 

0<t<* 



d? 



fk,i (t) 



3 



On the other hand if h is a function polyharmonic of order TV in the ball 
Br, then we have as in ([3]) the expansion 



ft (*) = (O»- fcy M(0). ( 5 ) 

k=0 1=1 

and it is well known (see Sobolev [5D], or [13]) that the coefficients hk,i{-) 
are polynomials of degree N — 1. Thus we may put into correspondence the 
functions fk,e and the polynomials hkj, which is the core of the polyharmonic 
interpolation. 

The polyharmonic interpolation problem is now very natural to formulate: 
Assume that for every fixed (k,£) with k = 0, 1, ... and £ = 1,2, dk we have 
interpolation points which we assume to be pairwise different: 

< r k ,t.\ < r k ,e,2 < •■■ < r k . LN < R. 

Then for every (k, £) we find the polynomials hk,i of degree < N — 1 from the 
one-dimensional interpolation problems 

hkt (r 2 kAj ) = he (r 2 kAj ) for j = 1, 2, N. (6) 

Now the main question is: For which distribution of the points \Xk,i,j} and 
for which functions / is the series in ([S]) convergent? If we have convergence 
then we will call the function h a polyharmonic interpolant of order TV. Our 
first result says that for every distribution of the points {vk/,j} and for a wide 
class of C°° functions / we have convergence. Indeed, we have the following 
amazing result. 

To make our result more transparent we will introduce the following semi- 
norms denoted by || • | jv"> wn i cn are motivated by 



|jy := lim sup 

M o<t<i? 



fd 



N 



1 



Nl \dt N 



fk,e (t) 



(7) 



We see that 



sup \\f\ N = M 

N>0 



where M is the constant in (|3|). 

Theorem 2 Let the function f be C°° in a neighborhood of the closed ball Br 
and the interpolation knots {rk,e,j} k i j satisfy 

< r MJ < R, 

where k = 0,1,2,..., £ = 1,2,..., dk and j = 1,2,...,N. If the seminorm \\f\ N 
satisfies 

R\\f\ N <l> ( 8 ) 



4 



then there exists a unique polyharmonic interpolation function h (r0) of order 
N which belongs to L2 for every r < R, and h belongs to L2 {Br) . 

Assuming Q), the error of interpolation is given by 

\\f(r9)-h(r6)]\ L2{sn - 1) <CR 2N \\f\% +1 . 

We see that in a certain sense the above Theorem [2] presents a complete 
analog to the one-dimensional interpolation since we may take arbitrary knots 
of interpolation Tk,e.,j- However we see that condition ((5J) is a restriction on 
the arbitrariness of the data / and this is the price which we have to pay 
for the infinite-dimensionality of the problem. Does this restriction imply a 
specialization in the one-dimensional case? The answer is "no". Indeed, since 
the one-dimensional polyharmonic functions of order N are just polynomials of 
degree < 2 N — 1 we see that condition (jHJ) is trivially fulfilled due to lim in ([7]) . 

There is still another way to consider the one-dimensional case embedded 
into the multivariate case, namely, when in the sums (J3]) and ([5]) only the term 
for k = is non-zero. Then / (x) = /o,i (r 2 ) and h (x) — ho,i (r 2 ) where /io,i (0 
is a polynomial of degree < N — 1. Indeed, in the one-dimensional a C°° function 
/ is identical with the univariate analytic function /o,i (•) in the expansion (|2|), 
and the knots are ^o.ij with 1 < j < N. We see that in this case restriction 
([5]) is always satisfied, i.e. Theorem [5] extends the one-dimensional theory in a 
natural way. 

If we change the point of view, and consider / to be fixed, then we have to 
choose a radius R small enough to fulfill (JSj) . 

As a second result we consider the special case of the knots which are lying 
on TV concentric spheres in R™, i.e. when the knots {rk.e.j} k t j satisfy 

r k ,i, 3 =rj for j = 0, 1, 2, iV - 1 

for all indices {k, £). Assume that / is a function analytic in a neighborhood of 
Br, and that the polyharmonic function h is an interpolant of /, i.e. satisfies 
([5]). From the expansions in spherical harmonics ([3]), ([5]) for every fixed r, and 
for j = 0, 1, 2, iV — 1, we see that the interpolation problem (JSJ) is equivalent 
to the following polyharmonic interpolation problem on concentric spheres 

h (rjO) = f ( rj 6) for e S"- 1 . (9) 

Let us recall the following result from Theorem XI. 3], 

Proposition 3 Let ip be a function defined and continuous on the unit sphere. 
A necessary and sufficient condition for the analyticity of tp is that in the rep- 
resentation 

00 d k 

k=a £=1 

the coefficients pu,i have exponential decay, i.e. there exist two constants K and 
n > 0, such that 

WkA < Ke^ k for every k = 0, 1, 2, t = 1,2, d k . (10) 
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Let us put 

tp* (9) := / (rjff) forOeS 71 - 1 . 

From the estimate (fTU|) we see that for all j = 0, 1, 2, N— 1 we have a number 
rjj > such that 

^{9)\<Ke-^ k . (11) 

Now we have again the question of convergence of the series ([5]) and it is 
solved by the second main result of our paper: 

Theorem 4 Let the numbers Tj with < T\ < r 2 < ... < rjy < R be given, and 
for the parameters rjj of the analytic functions tp 1 defined in ill]) the inequality 

R ■ max ( — - J < 1 (12) 

be satisfied. Then the polyharmonic function of order of polyharmonicity N 
satisfying the interpolation problem has an Li-convergent series in the ball 
B R . " 

Finally, let us remark that the polyharmonic interpolation problem |6]) may 
be considered as embedded in a more general scheme of interpolation theory [5] 
in the following way: Let us introduce the functionals 



L kAj if) = / / (tvj0) F M (0) d6. 

Then the polyharmonic interpolation problem ^ may be reformulated as the 
problem of finding the polyharmonic function h satisfying the infinite number 
of equations 

Lk.ej (h) = L kii j (/) for all k, £, j. 

On the other hand we have a nice demonstration of the polyharmonic paradigm 
|15j in the present situation. As we said in the introduction, the expectation 
that the knots Xi,x%, ...,xjsr in the one-dimensional interpolation theory will 
be replaced by closed surfaces jj, j = 1,2, ...,N in R ra in the polyharmonic 
interpolation has failed. Let us consider the sets 

Tj :={((kJ), PkAj ) :fc = 0,l,2,...; £ = 1,2,. ..,4} 

with pk,e,i < Pk,i,2 < ■■• < Pk,i,N ■ They may be considered as a multivariate 
generalization of the knots x\ < xi < ... < xn in the univariate case where Xj 
is replaced by Tj. For a better understanding of the role of the sets Tj let us 
make analogy with R n where the boundary 3D of a star shaped domain D in 
R ra (centered at the origin 0) can be written in spherical coordinates as 

3D = {(9,p g ) : for all 9 e S™" 1 } 

for some function pg > defined on the sphere §™ _1 . The results of the present 
paper show that the knots of interpolation x\ < X2 < ... < xn in the one- 
dimensional interpolation theory have been replaced by the sequence of mono- 
toncly increasing "sphere-like" sets Ti, T2, Tn- 
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2 Proof of Theorem [2t Poly harmonic Interpolant 
for General Knots 



Here we provide the proof of Theorem [5] 

Proof. By the definition of h in {5| and © hk,t are polynomials of degree 
< N — 1 and we may apply the classical results about the remainder of the 
interpolation, hence 



Rn-i (x) = 



(13) 



see [5] or [T71 (3.2.10)], and we obtain the formal series 



r 2 



here as usually w kt i (r 2 ) = Iljti (j ~ ' k,i,j 

By the definition of \\f\ N it follows by a standard argument that the L2 
norm of the above is estimated by 



\\f(r6)-h{ri 



Il 2 (S"- 1 ) 



dk 



\f(r6)-h(r9)\ 2 d9 



k+N+l 
N 



k=0 1=1 

AN 11 r\2N+2 ■ 



< EE K«« (O [I/ft 

: EEk fe n/i 



(14) 
(15) 



fc=0 £=1 



The convergence of the last series follows from the assumption ([8]). Hence 
follows the L 2 -convergence of the series for the polyharmonic function h. Also 
the estimate for the error of interpolation follows directly. ■ 



3 Proof of Theorem K Polyharmonic Interpola- 
tion on TV Concentric Spheres 

Next we prove Theorem 01 

Proof. Let us write the expansion of ip J in spherical harmonics 

00 d k 

^'(0) = EE<^w- as) 

fe=0 1=1 

Since the polyharmonic function h interpolating (p J on the sphere of radius 
Tj has the form 

h(x) = J2 Y ki (0) r k h M (r 2 ), 
k,t 
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where hki are polynomials of degree N — 1, we see that for all k = 0,1,2,. 
and I = 1, 2, <4, and for all j — 0, 1, 2, TV — 1 we need to have 



hence, 



h k i (rf ) 



We have to prove that the series for h is Z/2-convergent i.e. 



E /VMr 2 )| 2 dr< 



First we will find estimates for all h^.i- We have the explicit representation 
for the polynomials hkt in the form given in Krylov [T71 p. 42] and Davis [5J p. 
33] , where we put Xj = r\ ^ ■ . Let us put for the Lagrange fundamental functions 



k,l j 2\ 



Then we have 



(r 2 - x ) ■ • ■ (r 2 - Xj_i) (r 2 - x J+1 ) ■ ■ ■ (r 2 - £;y-i) 
(Xj - X ) ■ ■ ■ (Xj - Xj-i) (xj - x j+ i) ■ ■ ■ (Xj - X N -i) 



N-l 



hu (r 2 ) = E -Y (- 2 ) 



(17) 



3=0 



Bearing in mind that < r\ < ■ ■ ■ < r/v-i, we obtain (with the same K for all 
j's and rjj's ), the following estimate 



AT-l 

\hu (r 2 )|<^Ehf ( 



'r , 



7*™ 

j=0 J 



here K\ > is a suitable constant and 5 := miiij=i,2,...,N-i (\xj — Xj-i\) . Hence 
we obtain the estimate 



EE 

k=0 1=1 



oo d k .ft 



'*<EE 



„2A' 



" l „9 iv x ■* e 



^E 



7* • 

J=0 J 



< 



< 



A" 



2 oo d fc 



2 ^>27V 



JV 



^ 2W ) EE E 

fc=0 t— 

oo 

E fcn 



if! 

3=0 3 



N 



R 



dr 



2k+l 



k=0 



3=0 3 



2k + 1 

2 R 2k+1 
2k +1 



(18) 
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To obtain the last inequality we have used the estimate 



d k < Ck n ~ 2 



for some constant C > which follows from (JTJ) . Putting 

M = max ' 



3 \ r 



we obtain the estimate 



N e -mk 



Y,—i-<C(N + l)M k . 



r 

3=0 3 



The convergence of the series in (|18|l follows from the assumption 



R ■ max ( I < 1 . 



3 \ r 



Remark 5 If 



i?max ) > 1, 

3 \ Tn 



then in general one may not expect that the series representing the polyharmonic 
interpolant h will be convergent. This will be shown by the following example. 



Example. We assume that for all j we have 

— =C, 



(19) 



so that C = maxj \ ■ From ([IT)) it follows that 



N—l 

hu (r 2 ) = C k £ 4* CO - 

3=0 



and hence 



r k h u {r 2 )\ dr = £ 2k / r 2k 
o 



R 



N 

Ek.i ( 
W 3 ( 

3=0 



kJ I 2\ 

'r 



dr. 



According to the basic properties of the Lagrange coefficients (see e.g. [TTJ p. 
42-43]) 

AT-l 
3=0 
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so we get 



Jo Jo 
Finally, for a suitable constant C\ > the inequality 




r k h M (r 2 )\ 2 dr = C 2k [ r 2k dr = C 



■2k 



2k + r 




■2k 



R 2k+1 
2k + 1 



holds true and the last series is divergent since Ci? > 1. The proof is finished 
using assumption (fT§|) . 

If assumption p^|) does not hold then we can see by standard asymptotics 
arguments that for large k we will have 



for a suitable Ci > 0. This proves the divergence of the series. ■ 
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